Floquet topological quantum phase transitions in the transverse Wen-plaquette model 
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Our aim in this work is to study the nonequihbrium behavior of the topological quantum phase 
transition in the transverse Wen-plaquette model. We show that under the effect of a nonadiabatic 
driving the system exhibits a new topological phase and a rich phase diagram. We define generalized 
topological order parameters by considering cycle- averaged expectation values of string operators in 
a Floquet state. 
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I. INTRODUCTION 

Nature admits transitions between different states of 
matter driven by thermal fluctuations, which are referred 
to as classical phase transitions^^. At zero tempera- 
ture, as is expected, the thermal fluctuations do not play 
any role, however, the quantum fluctuations can entirely 
drive a change of phase in the system, which corresponds 
to a quantum phase transition (QPT). fn contrast to clas- 
sical phase transitions, a QPT is characterized by a dra- 
matic change of the ground state propertied. 

Symmetry breaking is a paradigm of condensed matter 
physics in which the states of matter are characterized as 
symmetry broken phases^ ^. Rather recently, novel states 
of matter have been found which cannot be classified 
inside the paradigm of Landau symmetry breakin^^^. 
These states of matter are characterized by an intrin- 
sic robustness against the effects of an environment as a 
consequence of the topological properties of the ground 
state. Examples of such topological states include the 
quantum Hall effect^ ^, topological insulators and topo- 
logical superconductors^^ This manifestation of crit- 
icality in nature is characterized by transitions between 
different topological quantum numbers, and corresponds 
to a topological quantum phase transition (TQPT). fn 
the integer quantum Hall effect, the different topolog- 
ical phases correspond to the Hall Plateaus observed 
experimentalljCl. A transition between states with dif- 
ferent Hall conductance can occur without symmetry 
breaking, fmportantly, further investigations in spin sys- 
tems with topological order will have a huge potential 
for applications in quantum information as topologically- 
protected qubits"''. 

Recent theoretical studies predicted the existence of 
topological insulators in solid state systemJ^^Jti^, with 
a subsequent ex perim ental observation in semiconductor 
heterostructureg^^^. Further investigations describe the 
emergence of Majorana modes at the edges of quantum 
wires^'' However, the experimental o bserva tion of the 
Majorana resonance is still under debat^^^E^l 

Despite of the striking results on TQPT in solid state 
systems, there is an increasing interest on the realization 
of topological states in other experimental setups. In 
particular, a realization of quantum magnets in a system 



of cold atoms placed in an optical lattice has been sug- 
gested recently^" •^■^ , to allow the implementation of the 
spin-1/2 Kitaev model on a hexagonal lattice. Further- 
more, topological protection can be realized by means 
of an array of superconducting nanocircuit^^H^ which is 
characterized by an Z2 topological order parametei!^. 

Previous investigations in QPTs inside the Landau 
symmetry-breaking paradigm, show that the external 
control induces effective interactions that lead to the ex- 
istence of new quantum phased^UM] ^ characteristic of 
the driven manybody systems is that the periodicity in 
time opens the possibility to create gapless excitations in 
a controlled wa y by t uning the parameters around quan- 
tum resonanced^maoi Furthermore, the influence of the 
external control on the topology of the system via geo- 
metric phases has also been explored"^^ 

The current experimental feasibilities have motivated 
the study of TPQTs under an external driving. Rather 
recently, the concept of topological charge of the Floquet 
Majorana fermions defined in terms of the Floquet op- 
erator has been introducecP^. Furthermore, under the 
effect of intense circularly polarized light, a gap can be 
open in the Dirac cone, leading to a photoinduced dc 
Hall current in graphen^^. Similarly, by using one- 
photon resonances, a trivial insulator can be driven into 
the top ologica l phase, originating a Floquet topological 
insulatoiIalEll. 

In contrast to the works previously mentioned, in this 
paper we show that the driving not only allows us to 
renormalize the parameters, such that the system en- 
ters into the topological phases, but, surprisingly, it in- 
duces a new topological phase that is absent in the un- 
driven system. In this work we study the Wen-plaquete 
model (WPM) in a monochromatic transverse field. The 
WPM has intriguing relations to other models studied 
in the literature. For example, it has been shown that 
the WPM can be exactly mapped to the toric code of 
Kitae\12£12II Furthermore, the toric code in a paralell 
magnetic field has low energy properties that resemble 
the two-dimensional transverse Ising model, which al- 
lows to study the influence of an external field in the 
TQPT'^^'^^. The effect of a transverse perturbation on 
the topological protection in the toric code has been ex- 
plored by using a mapping onto the Xu-Moore model. 
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which subsequently can be mapped onto the quantum 
compass model^. 

The paper is organized as follows: In Sec. |TT]we dis- 
cuss the general aspects of the transverse Wen-plaquette 
model (TWPM). In Sec. Ill we describe the physics of 



the system by means of the rotating wave approxima- 
tion (RWA) and discuss signatures of criticality based 
on the description of the nonequilibrium Ising transition. 
In Sec. |IV| we define generalized "string"- like topological 
order parameters by considering cycle-averaged expecta- 
tion values of string operators in a Floquet state. Finally, 
a discussion of the results is presented in Sec. W\ 



II. GENERALITIES OF THE TRANSVERSE 
WEN-PLAQUETTE MODEL 

In this section we describe the driven Wen-plaquette 
model in a transverse field on a, N x N square lattice 



Hit) = -git)J2Aj-jJ2^' 



where 



+ 1^4,^ + 1 



(1) 



(2) 



is the plaquette operator, and g(t) — go + gi cos57t. In 
this paper, Xi j,Yi j, and Zi j denote the Pauli operators 
acting on the (z,j)-th site of the square lattice. Fur- 
thermore, we assume periodic boundary conditions for 
even N . Fig. fl] (a) depicts the geometry of the plaque- 
tte operator, plaquette is even (odd) if the relation 
{-l)i+3 = 1 = -1] is satisfiecPS, as depicted in 

Fig.[l](6). 

In the topological phase, the undriven TWPM accom- 
modates three types of nonlocal excitations above the 
ground state of the system: Z2 charges, Z2 vortices, and 
fermions. Correspondingly, to describe these quasipar- 
ticles, the string operators of types Ti, T2, and are 
define(P3. The Ti string operator for a Z2 charge is given 

by 



(3) 
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FIG. 1: The transverse Wen-plaquette model. In (a) we 
represent the geometry of plaquette operator at the site («, j). 
Figure (b) depicts the real space lattice and the corresponding 
even and odd (— ) plaquettes. (c) Shows the Ti string 
operator along the path C, and the T2 string operator along 
C . Figure (d) shows the dual lattice and represents the A'^ 
decoupled Ising chains. To take into account the periodic 
boundary conditions, in the dual lattice the colors represent 
the different Ising chains. 



A. The duality transformation 



In this section, following the methods of 
we show that the driven TWPM Hamiltonian Eq. Q 
can be mapped into a set of decoupled one-dimensional 
Ising chains via a duality transformation. 

To obtain the mapping, we observe that the operators 
and Xi^j satisfy the commutation relations^ 



[F:,,F:A^[x,^b.Kd] = Q- (5) 



Wc is the product of spin operators along a path C con- 
necting even plaquettes of neighboring hnks (see Fig. [l] 
(c) ). Correspondingly, the T2 string operator for a Z2 
vortex is given by 



Wv(C) 



(4) 



In this case, is the product of spin operators along a 
path C connecting odd plaquettes of neighboring links 
(see Fig. [T (c)). Finally, the T3 string operators are 
defined as bound states of Ti and T2 strings, they are 
charge- vortex composite objects, and their definition is 
given in Refs.l^i^. 



Now, let us define a representation of this algebra by 



lattice 



F^ ; ^ 



and r/. on the dual 



(6) 



We parametrize the sites of the dual lattice by means of 
[i* = {i — j + 1, j) following the notation of Ref.l^. 
Under this duality transformation, the Hamiltonian Eq. 
([T| can be written in terms of N decoupled 1-dimensional 
Ising chains of length jv l37 | 4i|42 ] ^.j^^g paper i* denotes 
the chain index and j* the sites along the correspond- 
ing Ising chain. In the dual representation, the TWPM 
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FIG. 2: Caricature of the spin textures in the transverse 
Wen-plaquette model on the torus, (a) Depicts the topologi- 
cally trivial spin texture, (b) depicts a topological state char- 
acteristic of the undriven model, and (c) depicts the dual spin 
configuration depicted in[l](d). Finally, (d) shows the driving- 
induced topological state (in the spin texture on the torus, the 
dots represent spins pointing to the positive z-direction). 



Hamiltonian Eq. ([T|) is written as 

N 

where 



i' = l 



(7) 



N N 

Hm = -g{t) rf.,^..rf.,^..+i - J J] rf.,^-. (8) 

is the Hamiltonian of the z*-th Ising chain. The geometry 
of the decoupled Ising chains is shown in Fig. [l] (d) and 
Fig. 1(c). 

Therefore, to understand the behavior of the TWPM 
under the effect of driving, we must discuss the quan- 
tum criticality in the driven Ising model Eq. ([8|. How- 
ever, o ne can take advantage of the Kramers- Wannier self 
dualitjj™ 



(9) 



X . z z 

r<j* 

and obtain a dual spin Hamiltonian 

AT AT 

Hr,{t) = -9{t) -jH (10) 

which corresponds to an Ising model in a time-periodic 
transverse field^^. 



B. TQPT from a spin-polarized state to topological 
ordered state in the undriven case 



In this section we review the more important asp ects o f 
the TQPT in the undriven TWPM following Refs.EEIl. 

The quantum Ising model exhibits a second-order QPT 
go = J from a paramagnetic phase into a ferromag- 
netic phaseA However, to describe the TQPT in the 
TWPM one should define global order parameters as ex- 
pectation values of nonlocal operators in the ground state 
of the system; 



$1 = lim 6<)Xu = lim ( 60 rf ,.t-^ 



(11) 



(12) 



lim (erf ictJ^+i) 



In the topologically-ordered phase, go < J, the correla- 
tion function of af exhibits long-range order. As a con- 

sequence $i = 0, and $2 (l - (t)^)* ^ 0- This 
implies a macroscopic proliferation of Ti and T2 closed 
strings. In the spin-polarized phase (topologically triv- 
ial phase), go > J, a, long-range order of the correlation 
function of r^^ indicates the condensation of open strings 

close to the TQPT, i.e., $1 - ^1 - (^^^^^ ^ and 
$2-0. 

By using the duality between the Ising QPT and the 
TQPT in the TWPM, one we concludes that for go < J, 
the paramagnetic phase of the r-spin model (ferromag- 
netic phase of the cr-spin model) corresponds to the 
topologically-ordered (closed-string condensation) phase. 
Correspondingly, for go > J the ferromagnetic phase 
of the T-spin model (paramagnetic phase of the tr-spin 
model) is dual to the spin-polarized phase of the TWPM 
(open string condensation)^^ Figure [2] (a) depicts 

the spin-polarized phase of the TWPM on the torus and 
figure [2] (b) the topological phase. 

From the last discussion, we conclude that the global 
order parameters $1 and $2 are dual''"'^ Therefore, to 
characterize the topological phase transition, only one 
of them is necessary. Consequently, in this paper we 
consider generalized order parameters to characterize the 
topological phases, and we will not discuss the order pa- 
rameter characterizing the spin-polarized phase. 
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III. FLOQUET TOPOLOGICAL QUANTUM 
PHASE TRANSITIONS 



As we mention in the last section, the topological 
phases of the undriven TWPM are closely related to 
the quantum phases of the Ising model. In the driven 
case, however, the Ising model Eq. (10 1 exhibits a novel 
phase^^. In this section we perform a description of the 
TWPM based on the RWA™. 



A. The rotating wave approximation and the 
effective Hamiltonian approach 

Let us perform a unitary transformation of Hamilto- 
nian Eq. ([ij into a rotating frame via the unitary oper- 
ator 



Umit) = exp Wmit)^Xij 



N 



Ui^ ^rn (^) 



N 



i.j I i*=l 

N 



exp I 16*™ (i) ^ crf._j. 



(13) 



where 6m{t) — m{Q,/A)t+^ sniVtt. In the rotating frame, 
the Hamiltonian is given by Hm{t) = U^{t)'HUm{t), 
where H = H{t) — i-^ is the Floquet Hamiltonian^. The 
explicit form of this operator is given by 



-{l-cos[40„(<)]}^ 

— sin[46',„(t)] XijFi+i j-Xi+i j+iZi 

— sin[46'm(t)] XijZi+i j-Xi+i j+ili 

(14) 



where Am = go — m(f2/4) is the detuning from the m- 
photon resonancd^. 



The Hamiltonian Eq. ( 14 1 can be expanded in Fourier 
series 

oo 

4nW- ^i"^exp(inm) 

n— — OO 
oo N 

= E E^-r^exp(inl]0, (15) 



n — — oo 2*^1 



where /i^™^ and are time-independent operators. 



We define our effective Hamiltonian as an average of 
the Hamiltonian in the rotating frame Hm(t) over a pe- 
riod T = 27r/J7 of the driving 



Jo J ,.^1 



(16) 



where /^^-f'^^ is the effective Hamiltonian given in Eq. 



( 20 1 . Under the conditions 



4gi 



(17) 



the fast oscillating terms in the Hamiltonian Eq. ( 14 1 can 
be safely neglected and the effective Hamiltonian /iq"'' 
governs the dynamics in the rotating frame. In the last 
equation, Ji{z) is the Zth-order Bessel function^^. 

The time-independent effective Hamiltonian in the real 
lattice reads 



(m) 



Xi^jYiJ^\,jXi^i^j^\Yi, 



Jy 



E 



i+1 



Xi,jZi+i,jXi+ij-\-iZi^jj^i 



(18) 



where the parameters ji"^ = |[1 + (-l)™Jm(^)] and 
Jy™'' = f [1 ^ (^l)™i/m(%^)] are effective anisotropics in 
the rotating fram^2^ 

Besides the plaquette operator defined in Eq. 
([2]), the effective Hamiltonian Eq. (18) also contains a 
driving-induced plaquette operator 



^i.jZi^ijXi^ij^iZij^i. 



(19) 



This new term is absent in the original undriven TWPM 
and corresponds to an effectivejnteraction originated un- 
der nonequilibrium conditions^Sl 



B. AC-driven Quantum Phase transitions in the 
Ising model 

Our aim in this section is to describe the phase diagram 
of the driven TWPM Hamiltonian Eq. ([ij in terms of 
the quantum pha ses of the driven Ising model Eq. (10 1 
studied in Ref!^^. 

The criticality in the rotating frame is described by 
means of the effective Hamiltonian 



AT 



N 



f (cr,m) . \ ^ X T(m) \ ^ z 



N 



_ j{7n) V y 

j- = l 



(20) 



which is derived explicitly in Appendix \K\ 

Similarly to the Ising modei^, the Hamiltonian Eq. 



( 20 1 can be written in terms of Jordan- Wigner fermionic 



operators Cfc and c\. as follows 



?{cT.m) 



fc>0 



+ E(-l)"^- (^) h^^J-k + ^^kCk), (21) 
fe>0 ^ ^ 

where a;^ = 2 J cos k and /j, = 2 J sin 

The Hamiltonian Eq. (211 can be written as a Hamil- 
tonian for free fermions 



f (cr,m) 



k>0 

with dispersion relation 



(22) 



Efc.m = 2W (A,„ - J COS ky 



4.91 
f7 



sin fc 



(23) 

Therefore, when the gap between the positive and neg- 
ative energies closes, the effective Hamiltonian exhibits 
signatures of criticality in the rotating frame. Based on 
the results of Ref.^^, we find that the effective Hamilto- 
nian Eq. ( 20 1 exhibits effective ac-driven Ising-like QPTs 
along the critical lines |A,„| — J . Via the duality trans- 
formations Eqs. Q and the Ising-like QPT in the 
Ising model corresponds to a transition between the triv- 
ial and the topological phases in the driven TWPM. 

The Hamiltonian Eq. (20) also exhibit ac-driven 
anisotropic QPTs when the effective anisotropy 7'™) = 
Jrn (^) vanishes, under the condition |A,„| < J. The 
anisotropic transitions occur between two ferromagnetic 



phases. For ji™'' > Jy'"'' , the system exhibits long-range 
magnetic order along z direction (ferromagnetic phase 



FMZ), while for for J\ 



('») 



< J, 



(m) 



the dual Ising model 



Eq. ( pO| ) is ferromagnetically ordered along the y di- 
rection(ferromagnetic phase ¥MY^^. Interestingly, the 
anisotropic transition in the dual Ising model corresponds 
to an ac-driven anisotropic transition between topologi- 
cal phases of the TWPM. Figure |3] depicts the phase 
diagram for the ac-driven TQPT in TWPM for m = 0. 
The white zones in the phase diagram correspond to the 
topologically trivial phase (effective paramagnetic phase 
of the dual Ising model Eq. ( 10 1) and are defined by the 
inequality J < |Am| < lAmaTTfor m = 0, where A^ax 
denotes the maximum detuning for which the RWA is still 
valid. In Fig. E] (a) we plot the the effective asymmetries 



ji™^ and Jy"'' as a function of the driving amplitude. 

In the undriven TWPM the topological phases are 
characterized by nonlocal order parameters. However, 
as one can see after the last discussion, there is necessary 
to quantify the entities arising under the effect of driv- 
ing. In the next section we define generalized topological 
order parameters. 
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FIG. 3: Quantum phase diagram of the Wen-plaquette model 
in a time-dependent transverse field g{t) — go+gi cos fit. The 
figure depicts the quasienergy dispersion ±ek,m of the dual 
Ising model for parameters at the Ising-hke critical lines, in 
the ferromagnetic phases, and along an anisotropic critical 
line 7*™' — 0. For this plot we consider J/Q = 0.01 and 
m = 0. 



IV. CYCLE-AVERAGED TOPOLOGICAL 
ORDER PARAMETERS 

We now calculate global cycle-averaged expectation 
values of nonlocal operators. In particular, we focus on 
the expectation values in the negative-quasienergy Flo- 
quet state of the driven TWPM 



N 



M'(-) t 



where 



^ k,m - ^ 

k>0 ' k>0 

exp(-iiVi?(f)t) 

fc>o ■ 



i,(±) . 



(24) 



'^tlt) (25) 



is a Floquet state of the i*-th Ising chain. In the last ex- 



(±) 



pression, Em' — ± Jq ^£k,m is the total quasienergy 
and Sk.m is the quasienergy dispersion defined in Eq. 
(23 1. Furthermore, as we describe in Appendix [A] the 



Floquet fc-eigenmodes of the i*-th. Ising are obtained by 
applying an unitary transformation back into the labo- 
ratory frame \^''^l,t) = % ^.(t)|yl*2), with 



^fc.-(*)lxr._ 

xi^m) = cos(/)fc^„i|l_fc, Ifc) - sin(/)fc,,„|0_fc,Ofc), 
xiyii) = sin0A;,™|l_A;, Ife) +cos(/)fe,™|0_A;,Ofe), (26) 



satisfying /i^r^lxi^j,) = ±efc,m|xi*j,). We have used 
the basis of doubly occupied |l_/c,lfe) and unoccupied 
|0_fe, Ofc) states of ±fc fermions. Furthermore, we consider 
only the subspace with vanishing total momentuirP. 
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The Bogoliubov angle 
tion 

tan(2(/)fc,„) = 



is determined by the rela- 



2A„ 



(27) 



In the limit go ^ J, the negative-quasienergy Floquet 
state Eq. ( 24 1 corresponds to a topologically trivial spin 
texture in the absence of driving [see Fig. [2] (a)]. In 
the dual picture of the cr-spins, such a state is a tensor 
product of stationary states of the Hamiltonian Eq. ( 10 ) 
with all the spins polarized along the x axis 



exp 



(-iiv^4-)t) 



N 



\0-k,Ok). (28) 



i' = l k>0 



In the last expression 



9l 



,PE 



^90 J 



(.90 + jy 



(29) 



is the ground-state energy of the undriven Ising modeP, 
where E[z] is the complete elliptic integral of the second 
kincPS. 



A. Topological order from long-range magnetic 
order in the dual Ising model 

We focus here on the calculation of cycle-averaged 
global order parameters. The generalized closed-string 
order parameters are defined as products of plaquette 
operators along the diagonal string of the real lattice. 

By using the spin dualities Eq. ^ and Eq. ([9| we 
define 



$2 = 



lim — 

R^oo T 



^Ki dt 



1 1=1 



1 

lim - / (cr^ icrj 1+;^) dt, 



lim - 

i?— >-oo T Jq 



lim — 



1 1=1 



dt. 



(30) 



Our definition of $2 is a generalization to driven quan- 
tum syste ms of the global order parameters defined in 
equilibriuniSM]^ However, the existence of $3 has no 
analogue in equilibrium and a nonvanishing value of it 
would give rise to a topological configuration induced en- 
tirely by the external control. 

Interestingly, the cycle-averaged order parameters cor- 
respond to cycle-averaged correlation functions of the 
dual 1-dimensional driven Ising model Eq. (10 1. Such 



correlation functions characterize the long-range mag- 
netic order under the effect of external driving. One can 
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FIG. 4: (a) Depicts the effective anisotropics. The orange 
curve corresponds to J™ and the blue curve to J™, (b) Shows 
the nonlocal order parameters in the rotating frame. The or- 
ange curve denotes the order parameter C^z and the blue 
curve depicts Cyy. (c) Shows the cycle-averaged order pa- 
rameters. The green curve depicts $2, and the violet curve 
represents $3. We consider a G {y,z}, and the parameters 
J/Q. = 0.01 and m = 0. 



calculate the cycle-averaged nonequilibrium spin correla- 
tion functions in terms of the well-known results for the 
XY model^s* 51. 

For a G {y, z}, let us define nonlocal order parameters 
in the rotating frame Cqq = limj^_^oo Paa 

(R), where 



(31) 



corresponds to a spin correlation function calculated in 
the rotating frame (see Appendix [B| , where 



(-) 



fe>0 



A/c,m / ' 



(32) 



\xi ™> is given by Eq. IM\ 
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We can explicitly write (as we show in Appendix |B]) 
the topological order parameters as follows: 



$3^|[j(™)C,, + jMcj, 



(33) 



where ji"*-* , Jy™'' are effective anisotropics as they appear 



in Hamiltonian Eq. ( 20 1 



We find that $2 = $3 = in the topologically trivial 
phase |A„i| > J. For \A™\ < J the long-range correlation 
functions are given by 



Czz — 



^yy ~ 



2(|7<'"'|)^^^ 
1+|7(")| 



2(l7""'l)^/^ 



1/4 



^(™) > 
^(™) < 

^(™) > 

(A^)^ " - ^("0 < 



1/4 



(34) 



where 7'™-' = (— l)™j7m (^) is the anisotropy parame- 
ter. Fig. [4] (c) depicts the nonlocal order parameters in 
the rotating frame as a function of the driving amplitude. 

As a limiting case, our description recovers the physics 
of the undriven model {gi — 0) in the particular case 



771 = 0, where Aq = <?o, 7^°^ = 1, Jz = J, and 4"^ = 0. 
Therefore, we obtain that the nonlocal order parameters 
satisfy $2 — $2 and $3 — 0. 

We plot the dependence of the order parameters along 
the ladder of topological phases (ferromagnetic phases 
FMZ and FMy from the dual Ising model depicted in 
Figj3| in Fig. |4](c). Our results show that even for arbi- 
trarily weak driving amplitude gi, it is possible to drive 
the system from the trivial phase into the topological 
phase. Surprisingly, in the case m > 0, the static parame- 
ters are far from the critical point of the undriven system, 
i.e., gQ — J. For example, in the m > case, exact res- 
onance Am = implies to work in the regime where the 
undriven system is in the trivial phase g^ = ^ J. 

The strong driving regime allows to explore a new as- 
pect of our approach: by tuning the amplitude of the 
external control, it is possible to bring the system into a 
new topological phase. The phase diagram depicted in 
Fig. [3]shows rich patterns of multicritical points with no 
analogue in equilibrium when |A,„ — J\ and --yf™) = 0. 



V. CONCLUSIONS 

We have discussed the nonequilibrium TQPT in the 
Wen-plaquette model in a time-dependent transverse 
field. Similarly to the undriven model, a set of highly 
nonlocal spin-duality transformations allows a descrip- 
tion of the driven TWPM in terms of the Ising model 
physics. The topological character of the TWPM re- 
quires a description of the TQPT in terms of nonlocal 



order parameters, which correspond to long-range corre- 
lation functions for a driven Ising model in the dual pic- 
ture. We have introduced generalized "string"-like topo- 
logical order parameters by considering cycle-averaged 
expectation values of string operators in a Floquet state. 

In the case of con ventional nonequilibrium QPTs, pre- 
vious workJ^^H^HSlHini show that even under the effect of 
driving, the quantum states of matter correspond to sym- 
metry broken phases in the thermodynamic limit, i.e., 
the superradiant state does not conserve parity symme- 
try as in Ref.l^. In contrast, a TQPT corresponds to 
a change of phase without symmetry breaking. Similar 
to Refs.'^'^ , our approach reveals a possibility to induce 
topological configurations of the system by driving the 
topologically trivial phase. We show, however, that the 
monochromatic driving not only renormalizes the critical 
point, but generates an additional topological phase. 

Our methodology is potentially interesting in the con- 
text of quantum simulation with cold atoms. A promis- 
ing experimental realization of the driven Wen-plaquette 
model could be achievable by using a Rydberg atom 
quantum simulatoil^^!. In this setup, to simulate the sys- 
tem it is necessary to construct a quantum circuit con- 
sisting of nonlocal gates that encode the interactions of 
the effective model. 
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Appendix A: The rotating wave approximation and 
the effective Hamiltonian of the Ac-driven Ising 
model 



To obtain the effective Hamiltonian Eq. (20 1, we per 



form a description of th e system based on the rotating 
wave approximatiorpSlMI, Let us perform a unitary trans- 
formation of Hamiltonian Eq. (10 1 into a convenient ro- 



tating frame via the unitary operator 



Ui*.m{t) = exp 



^{t) > a: 



I exp 



fc>o 




-2z0„(<)(4cfc 



(Al) 

where 6',„(i) = m(ri/4)t -|- |y sinfJi, and are Jordan- 
Wigner fermions. 
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In the rotating frame the dynamics is governed by the and the Hamiltonian Hf, it) was defined in Eq. ( 10 1. The 
new Hamiltonian ifj'!^ m(0 = [^^i*,m(0]^^i*^j*,m(^)- The explicit form of the Hamiltonian in the rotating frame is 
operator Ui, = Hf. (t) - i ^ is the Floquet HamiltoniaiPl gi^^n by 



JV N N 

i?f.,„(t) = -A™ J2 - 2 + cos[40„(i)]} ^^'.3'^t','+i - 2 " E 



J ^ J ^ 



(A2) 



where A,„ — go — m(f2/4). By using the identity 

oo 

exp(izsinf2t) = ^ Ji{z) e'xp{ilQ,t) , (A3) 



chain Hamiltonian Eq. (10 1 



In particular, we are interested in the spin spatial cor- 
relations 



where J7i(z) is the Zth-order Bessel functiorP^, the Hamil- 
tonian Eq. ( A2 ) can be written in the form 



vco.{R,t) = (<^Ll.{t) (Bi) 



J2 hl:fe^p{^nnt). (A4) 



for a e {y, z}, where 



In analogy with the standard RWA of quantum optics, 
we obtain an approximate Hamiltonian to describe the 
mth resonance by neglecting all the terms in Hm (t) with 
oscillatory time-dependence: H[, ^(t) w h^"'™\ This 



W) = eM-^NE^-k)U,.,„Xt)\xV) (B2) 



approximation is valid as long as the conditions Eq. ([17|) is the negative-quasienergy Floquet eigenstate. 



hold. Finally, we obtain the time-independent effective 
Hamiltonian Eq. (|20l. 



In the last expression, 



(-) 



!q i|£fc,m is the to- 



tal quasienergy, and efe,m is the quasienergy dispersion 



Appendix B: Cycle-averaged correlation functions: 
Ac-driven long-range magnetic order 

In this section we perform the calculation of the spin 
correlation functions in nonequilibrium of the i*-th Ising 



defined in Eq. (23 1. Furthermore, the state \xm ') is de- 
fined in Eq. ([32 1 and the transformation Ui* rn{t) is given 
by Eq. dAl) . 



Without loss of generality, let us perform here the ex- 
plicit calculation for the particular case a = z: 



vMt) = (1 + cos[A9„rm (x^J xt^) + (1 - cos[49„M) (xt^ 

+ sin[40„,(i)] (xin^ \<i^Ii+b\ Xt^) - 

r 



M^o^{t)]{xi-^Ki<jli+R\xL^ 



(B3) 



After averaging over a cycle we get 



T 



(-) 



1 - 



451 



'14 1,1+-R 



(B4) 
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